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Aiming to describe frustrated quantum magnets with non-magnetic singlet ground states, we have extended 
the Rokhsar-Kivelson (RK) loop-expansion to derive a generalized Quantum Dimer Model containing only 
connected terms up to arbitrary order. For the square lattice frustrated Heisenberg antiferromagnet (J1-J2-J3 
model), an expansion up to 8th order shows that the leading correction to the original RK model comes from 
dimer moves on length-6 loops. This model free of the original sign problem is treated by advanced numerical 
techniques. The results suggest that a rotationally anisotropic plaquette phase |1| is the ground state of the 
Heisenberg model in the parameter region of maximum frustration. 



Over the last decades theoretical efforts have been devoted 
to study new quantum phases of bidimensional frustrated 
quantum magnets, motivated by the discovery of experimental 
antiferromagnets showing the absence of long-range magnetic 
ordering down to very low temperatures.- ^ * In such sytems 
a gap to magnetic excitation traditionally opens up while the 
spin-SU(2) symmetry remains unbroken. Two classes of "sin- 
glet" phases have been distinguished: Valence Bond Crystals 
(VBC) where some spatial symmetries are spontaneously bro- 
ken, and Spin Liquid (SL) for which all symmetries remain 
unbroken (e.g. the resonating valence bond (RVB) liquid.^). 

However, it is usually difficult to characterize the singlet 
phases in simple microscopic S = 1/2 models. For example, 
in the well-known Ji — J2 Heisenberg 5 = 1/2 antiferro- 
magnet on the square lattice, where frustration is controlled 
by the next-nearest interaction J2, no definitive answer has 
been given on the nature of the non-magnetic GS for maximal 
frustration at J2/J1 ~ 0.5. M. Mambrini et al. recently ad- 
dressed a work to this task^, studying the Ji — J2 — J3 model, 
containing an extra next-next-nearest neighbor J3 frustrating 
antiferromagnetic. 
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Interestingly, in this model the description in terms of Nearest 
Neighbor Valence Bond (NN VB) coverings^ is excellent in 
some extended region of parameter space, in particular around 
the point J2 = J3 1/4, with some significant extension 
along the line ( J2 + J3) = Ji/2. This model is therefore one 
simple canonical case where a truncation within the nearest- 
neighbor singlet configuration basis is legitimate and can be 
used as a simpler and convenient framework. 

In this paper, we have extended the Rokhsar-Kivelson 
(RK) loop-expansion to derive a generalized Quantum Dimer 
Model (QDM) acting in the space of hardcore NN dimer cov- 
erings of the lattice. We show that this expansion based on the 
hierarchy of the overlap matrix elements between the dimer 
coverings leads to an effective Hamiltonian that contains a 
sum of dimer moves, each involving only a single closed loop 
or loops at finite distances (connected term). In other words, 
all disconnected terms cancel out order by order We ap- 
ply this procedure to the Ji — J2 — J3 model and show that 



the leading contributions are of the form of a simple general- 
ized QDM on the square lattice which, in addition to original 
QDM'°, contains an additional loop-6 term which brings ki- 
netic competition in the system. The effective Hamiltonian 
then reads: 



eff 



(2) 



where the sums run respectivelly over all square or rectangu- 
lar plaquettes of the square lattice. Here we use the following 
convention: (i) White plaquettes denote kinetic (off-diagonal) 
operators that flip dimers around the thick contour (ii) Yel- 
low plaquettes stand for potential (diagonal) operators that 
leave configurations unchanged with a factor 1 if it is flip- 
pable around the thick contour and in the opposite case. In 
the following, the ommision of p indices in the diagrammatic 
notation is a shortcut notation for an implicit summation over 
all inequivalent plaquettes with a given shape. For example. 

Guided by a variational approach and combining numeri- 
cal techniques such as Exact Diagonalizations (ED) and Zero- 
temperature Green function Monte-Caiio (GFMC), we com- 
pute the phase diagram of this model. More specifically, we 
give evidence in favour of a Mixed Columnar-Plaquette phase 
first proposed in 1 1 ) and, since, evidenced in number of other 
contexts". Remarkably, this phase is found to be stable even 
in the presence of the loop-6 fluctuations. Hereafter, using the 
relation between the effective and microscopic models, we ar- 
gue in favor of the SU(2)-invariant version (i.e. applicable to a 
spin-1/2 model instead of a dimer model) of the above Mixed 
Columnar-Plaquette phase in the Ji — J2 — J3 microscopic 
model along the maximally frustrated line J2 + J3 = Ji/2 
where an approach restricting to the short-range VB basis has 
been justified previousIjEl. 

Derivation of the model. A systematic way to derive the 
generalized QDM hamiltonian ([T]) consists in (i) Projecting 
the Heisenberg model in the manifold formed by NN VB cov- 
erings of the square lattice, (ii) Perform the unitary transfor- 
mation that turns the non-orthogonal VB basis into the or- 
thogonal Quantum Dimer basis. The key ingredient of the 
calculation is the overlap matrix Oi^,^ = {^\ip) where \lp) 
and lip) are two NN VB states. Step (i) is equivalent to 
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TABLE I: Expansion of O, H and H up to order a^. Note that 
the effective hamiltonian contains only local (connected) processes. 
Some processes (marked as 0) does not appear in O or Ti, but are 
produced in Ti.'^^' by contractions of the generically non-commuting 
terms of the expansion (see appendix[B]for details). 



solving the generalized eigenvalue problem ^{{(p) — EO\(p), 
while the orthogonalization required in step (ii) is conve- 
niently achieved by defining Ti"''' = 0-1/271:0-1/2. 

Using a convention where all bond singlets are oriented 
from sites A to sites B according to the canonical biparti- 
tion of the square lattice, the overlap matrix can be written 
as O^^^ = Q,w-2ni(v» where N is the size of the system, 
ni{ip, ip) the number of loops in the ovelap diagram obtained 
by superimposing the two configurations, and a = l/\/2. 



On the other hand, {ip\Si.Sj\ip) = e((p\^p) with e = -3/4 
(resp. e = +3/4) if i and j belongs to the same loop of the 
overlap diagram but belong to distinct sublattices (resp. be- 
long to the same sublattice) and e = if i and j belongs to 
two distinct loops. Using a convenient scaling and shifting 
n (4/3)7i + JiN/2 of the Hamiltonian the matrix 
elements {(plHlip) can be expressed as Tii^,^ — h^^^O^^^ 
where h^p,^, only depends of the loops configuration. In par- 
ticular, this convention enforces h^p,^ = for all Lp. 

It is then possible to expand O and Ti in powers of a and 
compute = O'^^^TiO^^/^ accordingly as shown in ta- 
ble (uji up to . The expansion up to as well further tech- 
nical details of the calculation are given in the appendix |B] It 
is worth mentionning two peculiarities of this expansion : (i) 
contrary to several previous approaches'""'^"^' our expansion 
is not controlled by the length of the loops, but by the actual 
amplitudes of the overlap matrix elements that only depend 
on the overall number of loops in the overlap diagrams, (ii) 
all non-local and disconnected processes appearing in both Ti 
and O cancel in the expression of TY*^*. 

Let us discuss the results of this expansion. When trun- 
cated up to order a^, we recover the usual hamiltonian ob- 
tained in |10| with v/\t4\ = = 1/2. Note that such a 
drastic truncation appears a bit pathological in the sens that it 
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TABLE II: Quantum numbers of the eigenstates collapsing towards 
the same degenerate GS for each of the ordered phases considered in 
this paper. When applicable, we used the standard notations for the 
irreducible representations of the dv and C2v point groups, whose 
elements are defined w.r.t. a plaquette center. Definitions of the F, M 
and K points in the Brilloin zone are given in Fig.|2] ( x 2) denotes an 
additional first excited level (denoted by * in the text) in the {M, Ai) 
sector. The states with momenta Qs = (±27r/3,0), (0,±27r/3), 
(7r,±27r/3), (tt, ±27r/3), Q2 = (±7r/2,7r), (7r,±7r/2), Q3 = 
(±7r/2, ±7r/2) are even under reflection w.r.t. the momentum di- 
rections. The degeneracy of the pure columnar or plaquette (mixed) 
phases is 4 (8) and it is 12 for the trimerized phase. 



does not capture any aspect of the frustration of the original 
model : u/|t4| is indeed independant of J2/J1. In the per- 
spective of a justification of QDM model from the Heisenberg 
model, non trivial effects emerge from order a^. Fuithermore, 
considering the last column of table (jlji it is quite easy to see 
that, in the maximally frustrated region of the phase diagram 
{'h + Js, ^ Ji/'i) where the validity of the NN VB approach 
have been established*^, only the 3 processes retained in ^ are 
dominant. Importantly, we find that t4 > and te > which 
enable the use of efficient stochastic methods not applicable 
to the original frustrated spin model which suffers from the so 
called "minus sign" problem. 

Variational analysis: We now turn to the investigation of 
the effective Hamiltonian (|2|. We start with some discus- 
sion of the expected VBC phases shown on Fig. [T] Reg- 
ular columnar and plaquette phases have been introduced 
in the context of the frustrated Ji — J2 model and of the 
QDM'^. More recently, an anisotropic mixed columnar- 
plaquette phase has been introduced'. We consider here the 
possibility of such phases which interpolate between the sim- 
ple higher-symmetry VBC (such as columnar or plaquette). 
Because of the presence of loop-6 dimer moves, we also 
consider the possibility of a trimerization of the columns of 
dimers. We summarize the quantum numbers of the degener- 
ate GS of the various VBC in table|II] This will be used further 
in this paper to analyze the low-energy spectrum of Hamilto- 
nian (|2]l. 

Before showing the results of an extensive numerical anal- 
ysis, we first start with a simple variational analysis. Indeed, 
variational ansatze for the VBC phases of Fig.[T]can be easily 
constructed as tensor products of resonating plaquette states 
(see appendix |A]i and the knowledge of their relative stabil- 
ity provides a useful guide for the numerical search of VBC 
(but is also subject to some artifact of the variational method). 
For convenience, let us map the two-dimensional parameter 
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FIG. 1 : (Color online). VBC states considered in this work. Gen- 
eralized anisotropic VBC states labeled by CPa (mixed columnar- 
plaquette phases) and T (trimerized phase) interpolate between the 
most symmetric limits shown on the figure. 



space on a sphere by expressing the Hamiltonian parameters 
in terms of two Euler angles 6 and (j), as v = cos(0) sin(6'), 
t4 = cos(0) 008(6*) and tg = sin((^). The variational phase di- 
agram (Appendix[A| in the {9, (/))-plane contains three phases; 
(i) a RK region, (ii) the well-known completely isotropic 4- 
site plaquette phase and (iii) a large domain covered by the 
trimerized VBC with a 6-site unit cell interpolating between 
the columnar phase and the pure 6-site plaquette phase. For 
zero loop-6 kinetic term {(f) — 0), the phase diagram is the 
same than the one obtained in^(for 9 > — 0.257r). Once 
is turned on, kinetic fluctuations due to 6-loop kinetic terms 
suppress standard 4-site plaquette and RK phases in the vicin- 
ity of the RK point. It is remarkable that both the RK and 
the plaquette phase are rather robust under the additional 
kinetic term. However, the variational approach overesti- 
mate the stability of the RK phase (in fact limited to a single 
point with algebraic dimer correlations) and of the trimerized 
phase (the corresponding wavefunction has more 4-site flip- 
pable plaquettes than plaquette counterpart). In contrast, no 
mixed anisotropic VBC phase is found. A careful numerical 
approach is therefore necessary. 

Numerical results: We now move to ED of clusters up to 
8x8 sites. More especially, we compute the lowest en- 
ergy spectrum in each symmetry sectors, using both transla- 
tions and point group symmetries. Our analysis is based on 
the symmetry classification of the tour of states in the {9, 0) 
plane. In other words, each symmetry breaking VBC phase 
is characterized by a finite degeneracy of the GS with a set 
of well-defined quantum numbers (see table) separated by a 
gap from the continuum. On a finite size cluster, the de- 
generate GS is split but a close inspection of the low-energy 
spectrum can provide informations on the VBC phase (if the 
cluster is large enough). For — (i.e. cj) — 0), previous 
results (extrapolated to the thermodynamic limit) show that 
a phase transition between the columnar phase and a mixed 
columnar-plaquette phase (in fact the CPi phase of this pa- 
per) occurs around 9 ~ — O.OStt. Such a phase can obtained 
via an in phase spontaneous dimerization in the direction of 
the columns of dimers of the columnar phase or, equivalently, 
via a spontaneous rotation-symmetry breaking of the pure pla- 
quette phase. To simplify the discussion, we describe here two 
representative set of parameters that contain these two phases. 
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FIG. 2: (Color online). Typical ED low-energy spectra on a 8 x 8 
cluster (the GS energy is set to zero) for sl) 9 — — O.Stt and b) 
9 = tan"^ (0.5) as a function of 0. Levels of special symmetries 
(see text) are highlighted as colored symbols, from bottom to top: 
(r, ^i) (corresponding to the GS), [M, Ai), [F, Bi), {K, Ai) and 
(M, ^i)*. The arrow indicates the level crossing that marks the limit 
of the region where the later levels correspond to the lowest excita- 
tions, (c) Brillouin zone and its high symmetry points F — (0, 0), 
M = (tt.O) and a: = (tt.tt). 



9 = atan{0.5) (v/ti — 0.5) and 9 = — O.Stt, for which we 
have computed, as a function of (j), the spectrum of the ef- 
fective model by full ED. The spectra (defined w.rt. the re- 
spective GS energies) are displayed in Fig|2] Special symbols 
have been used to label five of the six low-energy levels (7 
states over 8) associated to the mixed CPi columnar-plaquette 
phase. We donnot consider the last one, the {K,Bi) level, 
which is believed to be more affected by finite size effects. 
The plots show wide intervals of (p where the above five levels 
are the true lowest eigen-energies, hence pointing towards the 
mixed phase as a possible GS. The level crossing at which the 
low energy spectrum becomes not anymore compatible with 
such a phase is indicated by an arrow. This level crossing can 
be used as a first crude estimator of the range of stability of the 
mixed phase. Surprisingly enough, the spectra at (? > is not 
drastically affected by a finite value of tg > 0), showing 
that the mixed phase is rather stable w.rt an extra loop-6 term, 
up to ~ O.Svr. This range of stability will be corroborated by 
our thermodynamic limit extrapolations of the order parame- 
ters (see below). In contrast, at 9 — —O.Stt and small (p, one 
can see that the very lowest levels of the spectrum (i.e. those 
really separated by a sizable gap from the rest) are compati- 
ble with the columnar phase whose symmetries are given in 
Table [III ^ narrow region of mixed columnar-plaquette phase 
might however exist at intermediate cj) values before the level 
crossing involving the {M, AI) state occurs. To finish this ED 
analysis, let us mention that, apart from the pure columnar and 
the mixed CPi columnar-plaquette phases, no region in pa- 
rameter space could be found where the low-energy spectrum 
is compatible with the other VBC phases described above. In 
particular, for large (relative) the spectrum becomes quite 
dense at low energies preventing any VBC phase identifica- 
tion 

In order to give a more quantitative determination of the re- 
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FIG. 3: (Color online). Extrapolations in the thermodynamic limit 
of the order parameters defined in the text characterizing the mixed 
CPi columnar-plaquette phase. Insets: finite size-scaling of both 
Af+(7r,7r) and M_(0,0) as a function of the inverse linear cluster 
size, using 36 sites (6 x 6), 64 sites (8 x 8), 100 sites (10 x 10), 144 
sites (12 X 12) and 196 sites (14 x 14) square clusters. The chosen 
value of 9 corresponds to the case of the Ji — J2 — J3 model studied 
in this work. 



gion of stability of the mixed CPi columnar-plaquette order, 
we have computed the related plaquette structure factors. 



{^o\Pp{-q)Pp{q)\^o) 



(3) 



where {q) is a diagonal operator with the same symmetry 
as the degenerate GS listed in Table. |ll] that we aim to tar- 
get, defined as Fourier transform of plaquette operators intro- 
duced in im. For positive and tg values, such quantities 
can be computed efficiently using Green Function's Quantum 
Monte Carlo (GFQMC). A Bragg-peak of I+{q) at momen- 
tum (tt, tt) (/f -point) and a divergence of /_ {q = 0) (F-point) 
reflect spontaneous translation and rotation symmetry break- 
ing of the mixed phase, respectively. Related order param- 
eters Mp{q) = y/Jp^q)/ L can be conveniently defined and 
results are displayed in Figj3] showing size-scalings of both 
Af+(7r, tt) and Af_ (0, 0) up to cluster size of 196 sites. These 
data correspond to the line v/t^ — 0.5, i.e. originating from 
the expansion of the microscopic model Ji — J2 — J3 under 
consideration here. Our results reveal that the Bragg peak at 
the ii' -point survives up to (/) < O.Stt, in good agreement with 
the ED criterium above. Interestingly, for increasing 0, the 
columnar order parameter /_ vanishes before the plaquette 
one: hence rotation symmetry is first recovered and a narrow 
region of pure plaquette order is stabilized between the mixed 
phase region and the more complicated (unknown) phase at 
larger t^. The extension of the GFQMC calculation to the 
9 > Q region which do not have ergodicity problems limita- 
tions, has enable to draw the phase diagram depicted in Fig|4] 
Concluding remarks: To finish, let us summarize our find- 
ings and their implications. First we have extended the QDM 
to the case with a finite tg amplitude for the loop-6 kinetic 
processes. Although we have also extended the search for new 
VBC phases, we found that the previously known columnar, 
plaquette and mixed (here called CPi) columnar-plaquette 




FIG. 4: (Color online). Phase diagram in the {9,(f)) plane obtained 
from numerical simulations. The previous knowledge of the = 
casJ^has been used, in particular to estimate the limit of the colum- 
nar phase as an approximate (black) straight-line. The boundary of 
the pure-plaquette iCP\ mixed) VBC phase obtained from the fi- 
nite size scaling of the associated order parameter is indicated by 
large circles (triangles). The (red) thick segment corresponds to 
the parameter region of the frustrated quantum antiferromagnet with 
J2 + J3 ~ Ji/2 according to the mapping described in the text. 
Crude ED estimates of the boundaries of the columnar phase and of 
the C'Pi phase are indicated by a thin dashed line and by small (blue) 
circles, respectively. The (light blue) region marked by a question- 
mark has not been identified. 



phases are stable when a moderate is added. This conclu- 
sion is first obtained by a close inspection of the low-energy 
spectrum of the model on finite clusters (ED). The quanti- 
tative determination of the extensions of the three previous 
phases is made possible by GFQMC simulations which do 
not suffer from the sign problem when tg > 0. Still, it has 
not been possible to characterize the GS in the whole param- 
eter space, in particular when ig dominates (accumulation of 
low-energy states) or when v < (i.e. 9 < 0) where the 
GFQMC encounters numeiical limitations (both cases being 
physically irrelevant anyway). Our present work on the gen- 
eralized QDM turns out to be very useful to make progress 
in the understanding of the frustrated S=l/2 quantum anti- 
ferromagnet on the square lattice. Indeed, it was previously 
argued that, in the vicinity of J2 ~ J3 ~ Ji/4 and along 
the J2 + J3 = Ji/2 line, a truncation in the NN SU(2)- 
singlet basis was legitimate. Using this argument and gener- 
alizing the RK expansion in terms of the magnitudes of the 
overlaps of the elements of the truncated (non-orthogonal) 
basis we have made a link between the microscopic model 
and some small region of parameter space of the general- 
ized QDM where, fortunately, a precise characterization of 
the phase can be made. We can therefore deduce that the frus- 
trated S=l/2 quantum antifeiTomagnet exhibits in the vicin- 
ity of J2 ~ J3 ~ Ji /4 the same type of lattice symmetry 
breaking as the mixed columnar-plaquette phase (CPi) of the 
QDM. In the language of the quantum antifeiTomagnet, it is a 
eight-fold degenerate SU(2)-symmetric phase with a 2 x 2 unit 
cell (like the plaquette phase) and rotation symmetry breaking 
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(like the columnar phase). While a previous investigation of 
the microscopic model on a 6 x 6 cluster has indeed provided 
evidence of plaquette correlation^, only the mapping to the 
effective model (which can be studied on much larger clusters) 
provides enough accuracy to show the spatially anisotropic 
nature of this spin-singlet plaquette VBC phase. 
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APPENDIX A: VARIATIONAL ANALYSIS 

The effective generalized Quantum Dimer Model on the 
square lattice originated from the microscopic frustrated 
Heisenberg Hamiltonian studied here can be first investigated 
by a variational method. Indeed, it is possible to construct 
variational ansatze for Valence Bond Crystal phases (see inset 
of Fig. |5]l as simple tensor products of resonating plaquette 
states, extending the number of possible phases arising in the 
standard QDM. Variational energies can then be computed an- 
alytically. Despite its simplicity, this approach reveals itself to 
be a useful guide for the numerical search of VBC (although 
artifacts due to its variational nature are expected) and can eas- 
ily incorporate the symmetry analysis provided in the paper 
In other words, there is a one-to-one correspondence between 
these variational wave functions (VWF) and the VBC states 
defined in the paper by the set of spontaneously broken sym- 
metries. 

We shall consider a set of five different VWF, (i) the well 
known \RK) one provided by Rokhsar and Kivelson'*' as an 
equal weight superposition of all dimer configurations, (ii) 
three VWF based on 4-site plaquette tensor products {\CPi), 
i = 1,2,3) and (iii) one with a 6-site unit cell interpolating be- 
tween the columnar phase and the pure 6-site plaquette phase 
(e.g. \T) in Fig.|5}. Excepted for \RK), all these VWF depend 
on a parameter a which allows a continuous interpolation be- 
tween pure singlet crystals and highly resonating VBC. To il- 
lustrate this construction, here is the expression of one of the 
anisotropic 4-site plaquette phases and the above mentioned 
6-site plaquette one: 
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FIG. 5: (a) Variational wave functions (associated to their respective 
VBC states) considered in this Appendix. Generalized anisotropic 
VWF labeled by \Aa) and j_B) interpolate between the most sym- 
metric limits shown on the figure, (b) Variational phase diagram as 
a function of 6 (x-axis) and (j) (y-axis). The thick red line shows 
the region of parameters of the effective model that maps to the 
Ji — J2 — Js model along the J2 + J3 ~ Ji/2 line for which 
the ground state is well described by singlet coverings. 



cos(q;) sin(a) for all the VWF, plus one in cos^(a) for the 6- 
site plaquette one. For the diagonal terms, both occupied and 
non-occupied plaquettes yield non-zero contributions to the 
expectation value of H. It is worth to emphasize that the \RK) 
WF requires a special analysis. Using well-defined Pfaffian 
techniques, one can compute analytically the probability of 
flipping a 4-site and a 6-site plaquette, which are respectively 
P4 = 1/4 and Pg = 0.0330112(1). Finally, these expectation 
values, as a function of the Hamiltonian parameters v{9,(j)), 
ti{0, (f) and ^5(^7 0)> are given by: 
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where the product is performed over the set of separate pla- 
quettes p or p' as suggested in Fig.js] 

In order to describe the stability of these phases w.rt. the 
parameters, the expectation value of the effective generalized 
QDM Hamiltonian is computed. The only off-diagonal terms 
of Ti. contributing to (r|7i|r), with |r) being one of those 
VWF, are plaquette flips on occupied plaquettes and the diag- 
onal potential term. This leads to contributions proportional to 



These energies are then minimized w.r.t. a and the varia- 
tional phase diagram displayed in Fig. |5] is obtained in the 
(6', 0)-plane. This phase diagram is discussed in the paper. 
The colors correspond to different values of the a parameter, 
a — tt/2 (blue) for the pure columnar phase, a = (purple) 
for the pure 6-site resonating plaquette phase and a = 7r/4 
(orange) for the isotropic 4-site plaquette phase. The RK re- 
gion has an arbitrary color. 




FIG. 6: (Color online). Overlap {tp\ij}) between two VB configura- 
tions If aAt}). Closed loops that appear in the overlap digram g are 
represented as yellow shades. 



FIG. 7: (Color online). The matrix element {(^|(4/3)Si.Sj|-0) ex- 
pressed as fij (y'l'i/')- Bond (i, j) is represented as a solid black line. 
Red and blue bonds represent \lp) and \tp). fij = +1 (respectively 
fij — — 1) if i and j belong to the same loop at even (respectively 
odd) distance along the loop, fij = if connects two sites on 
distinct loops (see Ref. il6i) . 



APPENDIX B: DERIVATION OF THE HAMILTONIAN 

This appendix is devoted to a technical description of the 
generaHzed QDM derivation scheme. 

Choice of a small parameter & disconnected processes. As 
briefly described in the article, the generalized QDM is ob- 
tained by developping = O^^^'^TLO^^/'^ according to 
the matrix element hierarchy of both O and Ti, in the VB ba- 
sis. Indeed, the amplitude of O^,^ (with <p and ?/; two VB 
configurations) only depends on the number of loops ni of the 
overlap diagram g = {tp, 4>) obtained by superimposing the 
two configurations : Og = a^-'^^-'ia) with a = l/\/2 and N 
the total number of sites (see figure[6|. 

The maximal number of loops is obtained for (p = i/j.ln this 
case g is just a collection of trivial loops with lentgh Ci 
2 and Og = 1. The next term is obtained by considering one 
non-trivial loop with length = 4 (^) while all remaining 
{N — 4) /2 loops are chosen trivial which leads to Og — a^. 
In the same spirit, one Ci = 6 loop i\ttt\) and (A^ — 6) /2 trivial 
loops is an a^-order process (see the fourth line of tables[l]and 
[lll| . Such a construction suggests that a quite natural driving 
parameter for the expansion is the length £/ of a unique non- 
trivial loop surrounded by {N — Ci)/2 trivial loops : such 
a process indeed appears at order a^'~'^. However the total 
length of loops is a constraint quantity {J2ieg ~ ^) 
other (non connected) terms indeed appear in O at the same 
order. For example, frtt formed by two disconnected squares 
also occurs in O with the amplitude a* despite the fact the 
non-trivial contour length is different from e.g. j±|:| . 

For this reason, in the derivation scheme presented here, we 
chose the overall number of loops as the driving parameter for 
the expansionof O rather than the commonly used length of 
the /oopJ^SEHIl The key difference lies in the presence of 
disconnected processes such as :ml : as we will see, they 
cancel at every order of the final effective hamiltonian, but are 
crucial in the calculation because they are responsible for the 
emergence of non trivial diagonal and off-diagonal connected 
processes. 

In the expansion of O, we use the following notation, 



p>0 



(Bl) 



where Up are combinations of uj^ process on graphs g : 



For a full list of lu^ up to order a^, see table 



III 



Processes 



Id 



n 



itjibitl: 



n 



-mm 
mm 



H 







2(J2- Ji)a" 



4(J2- Ji)a^ 



2 (-2Ji + 2J2 + Ja) a" 



6 (Ja - Ji) Q*^ 



2 (-3Ji + 3J2 + Ja) a'> 



2 (-3Ji + 3J2 + 2J3) 



2 (-3Ji + 3J2 + 2J3) 



8(J2 - Ji)q** 



2 (-4Ji + 4J2 + J3) 



4(-2Ji +2J2 + J3)q* 
4(-2Ji +2J2 + J3)a» 



4(-2Ji +2J2 + J3)a» 



2(-4Ji +4J2 + 3J3)q^ 



2(-4Ji +4J2 + 3J3)a^ 



2(-4Ji +4J2 + 3J3)a^ 



2(-4Ji +4J2 + 3J3)a^ 



(B2) 



TABLE III: Expansions of O and Ti. up to order a . 

Heisenberg hamiltonian expansion. The action of each 
term of the Heisenberg hamiltonian ([T]i of a VB state consist 
in a reconfiguration of dimers. Thus, ((^|7i|^) can be deduced 
form inspecting the topology of the overlap diagram {(p\4') as 
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Processes 




Analytic 
expression 


J2 = Ji/2 

J3 = 


J2 = Ji/4 
J3 = Ji/4 


J3 = Ji/2 
J2 = 




□ 


i-»/T r \ 4/'-ii 4\ 

2 (Ji - J2) a (1 + Q ) 


0.625 


0.46875 


0.3125 




□ 


-2 (Ji - J2) q2 (1 + Q^) 


-1.25 


-0.9375 


-0.625 


t 




-2 ((Ji - J2) (1 + fa^) - J3 (1 + 


-0.320312 


-0.304687 


-0.289062 


X 

4- 




2 (-J, + Jo + J-i) 


-0.125 


-0.125 


-0.125 


1 


II 


(-Ji + J2 + 2J3)a'' 


0. 


-0.03125 


-0.0625 




1 1 


( Ji - J2 - J3) q'' 


0.0625 


0.0625 


0.0625 


± 




(3Ji - 3J2 - J3) 


0.15625 


0.125 


0.09375 


I- 


■□ 


1 (-5Ji + 5J2 + J3) a* 


-0.140625 


-0.109375 


-0.078125 


I- 


■HI 


4(Ji - J2)q* 


0.25 


0.1875 


0.125 


i 

L 


J_L 


I (J, - Jo - J3) a* 


0.0234375 


0.0234375 


0.0234375 


P 




I (2Ji - 2J2 - 3J3)a* 


0.0078125 


0.0117187 


0.015625 






I (6J1 - 6J2 - 5J3)a* 


0.0546875 


0.0507812 


0.046875 


1 


ill II 


|(7Ji-7J2-3J3)a* 


0.0859375 


0.0703125 


0.0546875 




i 1 


(3Ji - 3J2 + J3) 


0.21875 


0.15625 


0.09375 




rrf 


-f (Ji- J2-2J3)a** 


0. 


-0.0117187 


-0.0234375 


IL 

n 


lH 


i(-5Ji + 5J2 + 6J3)a** 


-0.03125 


-0.0351562 


-0.0390625 






i(3Ji-3J2+4J3)a'* 


0.15625 


0.101562 


0.046875 



TABLE IV: Expansion of H up to order a*. 



recalled in figure |7] This allow H to be expanded in power of 
a (see Ref. 16^ similarly to O : 



E 

p>0 



(B3) 



Note that it is convenient here to rescale the hamiltonian 
by a factor 4/3. Furthermore, evaluating {ip\{4/3)Si.Sj\ip) 



generically involves an extensive number of lentgh-2 loops 
which only effect is to produce a trivial extensive contribution 
to the matrix element. This can be removed by scaling and 
shifting Ti to (4/3)7Y + JiN/2. Using this convention, the 
expansion of H contains only kinetic terms. For a full list of 
up to order a^, see the last column of table 



III 



Fusion rules & effective hamiltonian. The first step to com- 
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pute the effective hamiltonian Ti.''^ = O-^/^HO-^/^ is to 
derive the expression of O^^^^. To achieve this, we use the 
formal expansion : 



O 



^ rfi 



r(i + T) 



fc>0 " 



(l + r-fc)r(l + fc) 



(o-iy 



(B4) 



Powers of O and products with H generically involve sym- 
metric products of diagrams (see tabl e [Till ) that do not com- 
mute, e.g {\tL\, W.] } = \#\'^ m + \-f±±\ 1^ \±t\ . Evalutating 
these products requires establishing fusion rules that (i) gov- 
erns algebraic properties of diagrams and (ii) generate, order 
by order, new diagonal and off-diagonal processes. We give 
below, the minimal set of rule that are needed up to order . 
Order 4 fusion rule : 



Note that contractions of diagrams not only produce poten- 
tial (diagonal) terms, but also non-trivial assisted off-diagonal 
operators, such as e.g. that flip a plaquette only if a dimer 
is present next to it. Special process also appear such as 'mt 
(respectively i±[n:| ) which simultaneously flip two neighboring 
plaquettes with parallel (respectively perpendicular) dimers. 



In table IV we summurize the result of the expansion up to 
order a**. Interestingly enough, all disconnected terms vanish 
after simplifying the product O^^^^TiO^^^^. The demonstra- 
tion of this generic property is beyond the scope of the present 
paper and will be presented elsewher^S] At this level, let us 
remark that this absence of non-local terms in the generalized 
QDM Hamiltonian is physically satisfactory and is a strong 
indication of the internal consistency of the derivation scheme 
presented here. 



{gg, Eg} = + [tg + 2E§eS 



Order 6 fusion rules . 



\##\} = 2 ^ + 6 [tiw] + 2 + fm\ 
\\n\M} = 2\#\ + 2\m\ + \m\ 

Order 8 fusion rules : 
1 



{0, 0} = [g + 2 + 2 gg 



{Eg,Eg^} = 2 \tmt- 



1 



4 fliffl + 2 | t 1 t 1 t 1 - | + IjiHiffI + \-mt\ 

\mt\ + 2 \mm\ + \^^ + 6 \ m-vrn \ 
2 \-mmi\ + + \^m\ + ^ffl 



{Hi, M\} = gSg + gg + 2 



2 m. 



i j j j I ' 



{\#\, \m#\} = 2 + 8 \tmmt 
{H, gSg} = ggg + 2 \tm\ + 4 



+ 2 + 2 |fFFFf| + 2 

{§g, \m\} = \m\ + \M\ + 2 \mm \ 



+2\^+ m + m 

f\+2 



■2 m 



{\tt\-, :|±H:| } = 4 :n':rt + 4 'rtri. + 4 :i±i±n: + Jriiirt 
+ 2 .M u. + .cn-i 
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